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Abstract 

Production of doubly charged Higgs bosons via the s-channel process e~e~ — > H — > l~l~ 
at the future linear collider energies is studied by taking into account initial state radiation and 



O 
Q ' beamsstrahlung (ISR+BS), final state radiation (FSR) and detector smearing effects. The discovery 

bounds of lepton flavour conserving and violating couplings are obtained for doubly charged Higgs 

bosons. It is found that the future linear colliders with the center of mass energies \fs = 500 GeV 

:^J | and yfs = 3000 GeV will be able to probe the doubly charged Higgs bosons with diagonal couplings 



down to 10 and 10 3 , respectively. 



I. INTRODUCTION 

There is a consensus in the particle physics community that the next large scale project 
should be a linear collider. The International Linear Collider (ILC) |l| is an e + e~ linear 
collider with a center of mass energy of y/s = 500 GeV and the luminosity L « 10 34 cm _2 s _1 . 
Beyond this baseline machine, upgrades and options are envisaged whose weight, priority 
and realization will depend upon the results obtained at the Large Hadron Collider (LHC) 
and the baseline ILC. 

The physics results obtained in the first few years of running with ILC, together with the 
results from the LHC will then define the schedule for upgrades or other modes of operations 
(options). The expected shutdowns to install the upgrades or options will not take more 
than two years after an initial physics running time of at least four years, including the 
commissioning of the upgrades or options (l|. The ILC can operate for part of its running 
time in the electron mode where certain signals can be remarkably free from backgrounds. 
The e~e~ running mode can also offer the high polarizations (at least 80%) for both initial 
states. This mode is a particularly interesting feature of a high energy linear collider. 
The main feature of this option is that the e~e~ initial state is doubly charged and caries 
double lepton number. In this mode, the standard model (SM) activity can be highly 
suppressed, and clean signals from any non-standard physics can be searched. One of the 
many promising processes which can be studied in this mode is the resonant production of 
doubly charged scalars. Here we should note that s-channel resonance production of doubly 
charged Higgs bosons at e~e~ colliders dominates over their t-channel indirect productions 
at e + e~ colliders. Therefore, an e~e~ collider could be used as a factory for doubly charged 
Higgs bosons, where their masses, total decay widths and couplings could be measured 
precisely. An additional feature often associated with doubly charged Higgs boson is the 
possibility of lepton number violation. 

In many cases detailed measurements at a higher energy e + e~ collider would be needed 
to complement previous exploratory observations. We think that further progress in particle 
physics will necessitate clean experiments at multi-TeV energies which is possible at e + e~ 
Compact Linear Collider (CLIC) |2|. It is anticipated that the center of mass energy of the 
CLIC as high as 5 TeV might eventually prove feasible. Thus, the pair production of new 
particles could be detected in the e + e~ mode of operation up to the mass of \fsj1. 



The Standard Model (SM) gives a good description of the known fundamental particles, 
using SU(3)c x SU(2)l x U(1)y gauge group to describe their color and electroweak in- 
teractions. The SU{2)l x U(1)y electroweak gauge symmetry is broken to U(l) em by the 
Higgs mechanism, but a Higgs boson has yet to be observed |3j. This is one of the good 
reasons that other symmetry breaking mechanisms and extended Higgs sectors have not 
been excluded from the theoretical point of view. Some theories beyond the SM predict 
the existence of doubly charged Higgs bosons, including in left-right symmetric models (j|, 
Higgs triplet models 5j, and little Higgs models (a). In the compositeness models they can 
be considered as the doubly charged scalar bileptons p, |8j. The left-right symmetric (LR) 
electroweak theory based on the gauge symmetry group SU(2)l x SU{2)r x U(1)b-l was 
proposed to offer a dynamical solution to the parity violation of weak interactions. The 
presence of triplet representations of Higgs fields, i.e. H R triplet of SU(2) R and H L triplet 
of SU(2)l, provides a simple explanation to the lightness of neutrinos via the see-saw mech- 
anism l9j|, consistent with recent data on neutrino oscillations InJl. In addition, indication 

n 

for neutrino oscillations necessarily violate lepton flavor symmetry [11|. It has been empha- 
sized that the see-saw mechanism could lead to doubly charged Higgs bosons with masses 
accessible to current and future colliders. 

Direct searches for dou bly charged Higgs bosons have been performed by the OPAL 
|l2j, L3 |l3j and DELPHI |]J| Collaborations at LEP, and they have excluded H ±± bosons 
below masses of about 100 GeV, assuming exclusive H ±:k decay to a given dilepton channel. 
Recent searches by the CDF and DO Collaborations at the Fermilab Tevatron in the \x\x 
channel have excluded H^H^) below a mass of 136 (113) GeV |lj and 118 GeV (lj| at 
the 95% confidence level (C.L.), respectively. There are unknown parameters on which the 
obtained constraints depend: the mass of the doubly charged Higgs boson m H ±± and the 
coupling constants hij, where i, j = e,fi,r. The H ±± couplings hy to electrons and muons 
are experimentally constrained by the absence of H ±± production in e + e~ collisions and 
nonobservation of the decay \i — > 3e (h PP h PI , < 3.2 x 10~ 11 m H -±± GeV -2 ) and decay /i — > ej 
{hfijjien < 2 x I0~ 10 m^j±± GeV -2 ) (l3|. From the Bhabba scattering the relevant bound is 
given as h 2 ee < 6.0 x 10~ 6 m 2 H ±± GeV -2 |la| and from (g — 2) /t measurements it is /i 2 M < 
2.5 x 10 -5 77rL ++ GeV -2 [l9j. More stringent bounds exist from the muonium-antimuonium 
transition |2(j| in the form of the product of the couplings, h ee h^ < 2.0 x lO~ 7 m 2 H±± GeV -2 . 
Similar bounds are given for both lepton flavor conserving couplings h 2 ee < 3.5 x 10 -5 m^ 



GeV -2 |21fl and violating couplings h 2 , T ) < 1 x 10~ 6 m^__ GeV -2 [22] of doubly charged 
scalar bileptons (as doubly charged Higgs boson H ) from Bhabba scattering with LEP 
data. The experimental constraints on the couplings weaken with increasing H ±± mass. 

It is reasonable to assume that below the scale of 500 GeV the left-right symmetry is 

due to the strength of the left/right handed interactions and the form of CKM 

matrix of the right-handed interactions. It is entirely reasonable to suppose that a H in 

the mil — < 500 GeV mass range relevant for a y/s = 500 GeV e~e~ collider would already 

have been observed at the LHC, regardless of the magnitude of the coupling hu. 

In e~e~ collisions the production of doubly charged Higgs bosons via lepton number 
conserving and violating processes have been studied in (a,l2J|- A work for the usefulness 
of looking at associated monoenergetic photon as signals for a doubly charged scalars at a 
linear e~e~ collider with yfs = 1 TeV has been presented in [25]. Recently, doubly charged 
scalar bileptons (as doubly charged Higgs bosons) have been studied via e - e~ — > e~e~ , \x~ \T 
processes with a resonance scan of the collider energies up to 2 TeV [26]. 

In this paper, we study in detail the production of doubly charged Higgs bosons H£ R 
through s-channel in e~e~ collision including initial state radiation (ISR)+beamstrahlung 
(BS) and final state radiation (FSR). We also include the smearing effects of a detector on the 
four-momentum of final state leptons. We analyze the signature for identifying an isolated 
same sign, planar and pr balanced two-lepton events. We also calculate the discovery bounds 
of couplings and masses for doubly charged Higgs bosons with lepton flavour conserving and 
violating couplings. 

II. DOUBLY CHARGED LEPTONIC INTERACTIONS 

The doubly charged Higgs bosons have the following Yukawa couplings to the leptons: 

L Y = h R>ij lf R Ca 2 H R l jjR + h Ltij lf L Ca 2 H L l jtL + h.c. (1) 

where lJ L t R \ = {^i,h)L(R) are two-component left or right handed lepton fields with flavour 
index i; C is the charge conjugation matrix; a 2 is the Pauli matrix. Hi and H R are the 



SU(2)l and SU(2)r triplets, respectively. The Higgs triplets 



H^\" l if 2 " Z7R ir \ (2) 

acquire non-vanishing vacuum expectation values (vev) given by 

1 ( , 

^^vsUo. 1 (3) 

where the left-handed triplet vev vl is kept to be small due to its contribution to the p- 
parameter, p ~ (1 + 2v 2 L / k 2 )/{1 + Av 2 L / n 2 ) where k is related to the vev of bidoublet scalar 
field. The vev vr gives masses to Wr as well as doubly charged scalars H^ R . The strength 
of the interaction is scaled by the symmetric unknown constants hij which, in general, are 
not flavour diagonal allowing for lepton number violating interactions. Therefore there are 
six possible couplings: 
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Using the Lagrangian (1) with the assumption on the structure of the leptonic couplings, 
the leptonic decay widths for the doubly charged Higgs bosons are given by 



, __ = m H-\^ u2 



in 

1,3 



D4 ( 5 ) 



We calculate the leptonic decay widths for H as Th — = 0.597 GeV taking the mass 
m H — = 500 GeV and assuming universal flavour diagonal couplings hij = 0.15^-. Since 
we consider only leptonic decay channels of doubly charged Higgs bosons, calculated decay 
widths become narrower for smaller couplings and masses. Other decay channels including 
gauge bosons can be neglected assuming the vev of the neutral member of the triplet is 
small and taking into account the limits on the model's couplings and the mass spectrum 
of the Higgs bosons. The two doubly charged Higgs bosons H~[ R ~ can have different chiral 
couplings to leptons. Their masses are expected to be comparable with each other because 
they derive both from similar terms of scalar potential. Since their production processes are 
the same we concentrate on the production of H (H£~ or H R ). 
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Figure 1: Feynman diagrams for the process e e — > e e including doubly charged Higgs boson 
H exchange. 

III. PRODUCTION CROSS SECTIONS 

For the process e~e~ — > e~e~ including doubly charged Higgs bosons H^ R , as depicted 
in Fig. [H the corresponding differential cross section can be obtained as 
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where Y H — is the total decay width of H L R boson. The polarization factors F RRLL = 
(l±Pi)(l ± P 2 )/4 and F RLiLR = (l-P 1 P 2 )/2 with the polarizations P x and P 2 of the in- 
coming electron beams. The mixing angle dependent constants Cll = s^/cj,, C RR = 
(1 — 2s 2 u ) 2 /(2s w c w ) 2 and Cl R = C R l = (1 — s 2 u )/2c 2 u . The abbreviations c w and s w are 
used for cos^jy and sinOw, respectively. In Eq. @, g e is the electromagnetic coupling 
constant; s and t are the Mandelstam variables, and mz is the Z boson mass. 

Here we assume that electron beam is ~ 100% left or right polarized. This is a simpli- 
fication which is motivated by the very high polarization rate 80 — 90% achievable at the 
e - e~ options of the ILC and CLIC environment. The impact of an exact amount of 80% 
left polarization results in the mixture of partial LL, LR/RL and RR contributions (being 
LL 80% dominant) to the cross section given in (6). 

For the processes e~e~ — > yTpT and e~e~ — > t~t~ only flavour diagonal couplings 
contribute to the cross section. The resonance cross section for the signal processes e~e~ — > 



H — »■ I I (where I =\i or r ) can be obtained as 

87iT ee r M s 

° = m 2 H - [(s-m%-)* + m%__T%__] {V 

if the decay width of doubly charged Higgs boson is significantly larger than the beam energy 

spread (typically 10~ 2 a/s). It is a good approximation to assume monochromatic electron 

beams. However, in the narrow width approximation (r# — jvan — < 0.01), the cross section 

formula (7) becomes 

a = J^B ee B u T H -5{s - m 2 H ^) (8) 

m H — 

depending on the total decay width and branchings Bu of the doubly charged Higgs boson. 
Here, we consider only leptonic decay channels of the H and assume that the energy 
spectrum of the luminosity is broader than the narrow resonant width. It should be noted 
that the actual total cross section can be obtained by convoluting a(s) with the initial and 
final state radiation functions, and the smearing function. 

We have implemented the doubly charged Higgs boson interaction vertices into the Monte 
Carlo simulation program CompHEP |2j| for decay width and cross section calculations. 
Since the center of mass energy is not well defined due to initial state radiation, beam- 
strahlung and the energy spread of the particles in the beam, these effects modify the 
energy dependence of the resonance cross section significantly. If the luminosity spectrum is 
well measured the physical parameters can still be extracted from the resonance scan. The 
scan can be done measuring the cross section at different center of mass energies around the 
resonance point. 

In the case of the process e~e~ — > H — > /i~ \T the cross section with and without 
ISR+BS are calculated around the resonance point as shown in Fig. 01 and El The resonance 
peak shifts to the right (a higher energy) and reduce by a factor of about 5 due to ISR+BS 
effects. At the ILC based e~e~ collider with the center of mass energy \fs = 500 GeV we 
obtain the resonance cross sections: a = 2.17 x 10 3 pb, uisr = 1-03 x 10 3 pb and <Jisr+bs — 
4.45 x 10 2 pb; without radiation, with ISR and with ISR+BS effects, respectively. For 
CLIC based e~e~ collider with the center of mass energy y/s = 3000 GeV this effect become 
apparent from Fig. as we find the signal cross section a = 60.41 pb without radiation 
effects. In this case ISR+BS affects the cross sections significantly as <Jisr+bs — 1-84 pb. 
When calculating these effects we take into account the beam parameters of the ILC and 
CLIC based e~e~ colliders as shown in Table HI 
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Figure 2: The production cross section for the process e e — > H — ► /j, \x depending on the 
centre of mass energy of the e~e~ collider. Solid line shows the cross section without the ISR+BS, 



dashed line represent ISR+BS effect for the flavour diagonal couplings h, 
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Table I: The beam parameters of ILC and CLIC based e~e~ colliders used in the calculations of 
ISR and beamstrahlung. N is the number of particles in the electron bunch; cr x ,y,z are the average 
sizes of the electron bunches; T is the beamstrahlung parameter and _/V 7 is the average number of 
photons per electron. 
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Figure 3: The cross section a(e e — ► H — > fi \x ) depending on the center of mass energy 
around 3000 GeV. Solid line for the cross section without the ISR+BS, dashed line represent ISR+BS 



effect for h eetM1 = 0.1. 



In Fig. B]and [13 we present the differential cross sections (da/dpT) for the signal (h ee = 
0.1) and background through the process e~e~ — > e'e" at the center of mass energy y/s = 
500 GeV and y/s = 3000 GeV. Here we take the mass of H '" boson equal to the center of 
mass energy of the collider. We see from these figures that the signal can be differentiated 
from the background by applying the appropriate px cuts on the electrons. 

Since we mainly focus on the resonance production of doubly charged Higgs boson via the 
process e~e~ — > e~e~ we may neglect the SM/H interference terms around the resonance. 
Because, we obtained the best limits for the couplings in this case. We see that after the 
Pt cuts these terms contribute insignificantly. Furthermore, in simulation with standard 
PYTHIA we can't produce signal and background events together and their interference at 
the same time. Therefore, we will take into account signal and background separately and 
make the statistical analysis accordingly. 

In Table ITU and Mil we show the signal and background cross sections, for a more quan- 
titative outlook, depending on the applied pr cuts for y/s = 500 GeV and y/s = 3000 GeV, 
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Figure 4: Transverse momentum distribution of final state electrons for the signal (h ee = 0.1) and 
background processes at y^ = 500 GeV. 
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Figure 5: Transverse momentum distribution of final state electrons for the signal (h ee = 0.1) and 
background processes at y/s = 3000 GeV. 
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Table II: Cross sections for the signal and background with and without ISR+BS effects depending 
on the applied pt cuts on final state leptons at ILC energies. 



Pr(GeV) 


<7 5 (pb) 


&S,ISR+Bs{pb) 


<j B (pb) 


<7B,ISR+Bs{pb) 


10 


2.17 x 10 3 


4.45 x 10 2 


3.00 x 10 3 


2.99 x 10 3 


50 


2.13 x 10 3 


4.36 x 10 2 


1.27 x 10 2 


1.264 x 10 2 


100 


1.99 x 10 3 


4.07 x 10 2 


3.41 x 10 1 


3.35 x 10 1 


150 


1.74 x 10 3 


3.56 x 10 2 


1.52 x 10 1 


1.44 x 10 1 


200 


1.31 x 10 3 


2.67 x 10 2 


7.29 x 10° 


6.28 x 10° 



Table III: The signal and background cross sections depending on the applied pt cuts on final state 
lepton with and without ISR+BS effects at CLIC energies. 



PT (GeV) 


<?s (pb) 


CS,ISR+Bs(pb) 


o-fi(pb) 


&B,ISR+BS(pb) 


100 


6.05 x 10 1 


1.81 x 10° 


3.23 x 10 1 


3.02 x 10 1 


200 


5.95 x 10 1 


1.79 x 10° 


8.46 x 10° 


6.11 x 10° 


500 


5.70 x 10 1 


1.71 x 10° 


1.44 x 10° 


4.74 x 10" 1 


1000 


4.50 x 10 1 


1.35 x 10° 


3.46 x 10" 1 


4.94 x 10~ 2 



respectively. The background cross section decreases significantly with the higher px cuts. 
Some optimized values can be found by searching for signal to background ratio. We have 
used a relation for the determination of optimal pt cuts pr > run — /2 — Ap T where Apr can 
be optimized for every mass values to give better signal to background ratio and resolution 
as well. With a cut pr > 200 GeV on the final state electron we obtain the background 
cross section Ob = 6.28 pb and signal cross section as = 267.0 pb for m H — = y^ = 500 
GeV and h ee = 0.1. In the e~e~ collision with y/s = 3000, if we apply a pt > 1000 GeV we 
obtain a smaller cross section for background Ob = 4.94 x 10 -2 pb and signal as = 1-35 pb 
for mu — = V^ an d h ee = 0.1. 
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IV. SIMULATION 

To make our analysis realistic we have smeared the energy of final state particles by a 

Gaussian function whose half- width is guided by the resolution of the electromagnetic (EM) 

calorimeter [28] 

SE 10% 

ir-_ +0 - 01 (9) 

Energy resolution defines the detector's ability to produce narrow peaks in the spectrum. It 
is affected by detector geometry and electronic noise. Resolution is typically measured as 
the width at the full width at half maximum for the given energy. In all cases the smaller 
the number specified for energy resolution is best and gives the user a better probability of 
separating peaks for easier identification. The invariant mass resolution of the final state 
charged leptons can be given by 5mu/mu ~ 5E/\^2E. Here, we calculate the invariant mass 
resolution 5mu = 4 GeV for m« = 500 GeV. When reconstructing the invariant mass of two 
leptons in the final state we have incorporated the effects of initial (ISR) and final state ra- 
diation (FSR) using the PYTHIA program |2flj |. We use the program SMEAR |2Qj to smear 
four-momenta of final state charged leptons, taking into account detector effects and accep- 
tances. We have included the detector effects in our calculations for ILC and CLIC energies. 
For electrons in the final state, the magnitude of the output momentum is determined from 
the smeared calorimeter energy assuming the mass of the particle is that of an electron. 
A tracking device is used for the measurement of four-momenta of the muons in the final 
state. The tracking device is modelled in the following way: (i) the transverse momentum 
Pt and azimuth angle are smeared together, taking into account their correlation [31]; 
(ii) the z-component of the momentum is smeared by projecting the transverse momentum 
according to a Gaussian in l/p z . A detailed discussion of the properties of the detector 
model implemented in the present version of SMEAR can be found in J3i|. Charged leptons 
(electrons and muons) are reconstructed within the acceptance range in pseudo-rapidity 
\r]\ < 2. Momentum resolution of the tracker is given depending on the scattering angle, 
and assumed as 0.02% for 6 > 40 degrees. The efficiencies for charged particles are assumed 
to be 98%. 

We show the invariant mass distribution of the e~e~ pair for the above choice of parame- 
ters in Fig. H3 Here the distribution peaks correspond to the mass values of doubly charged 
Higgs bosons. Invariant mass of two electrons is smeared according to the detector effects 
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Figure 6: Invariant mass distribution of two electrons at yfs = 500 GeV for h ee = 0.1: (a) signal as 
a narrow peak, (b) signal with detector effects and a Gaussian fit. 

to be observable experimentally. 

The invariant mass distribution of two electrons in the final state at the ILC and CLIC 
energies are shown in Fig. and [HI Since we include the ISR+BS and smearing effects 
to the production cross sections there appears a coupling dependence moderately at very 
resonance point. In Fig. we present the signal with m H — =350, 400, 450, 500 GeV 
and the background (e - e~ — ► e~e~) with p Te > 100 GeV and pr e > 200 GeV. In Fig. |Hl 
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we present the signal with m H — =500, 1000, 2000, 3000 GeV and the background with 
pre > 200 GeV and pr e > 500 GeV. We see from these figures that the signal peaks over 
the background at resonance even for smaller couplings. Applying appropriate Pt cuts on 
the electron from the background process (e~e~ — ► e~e~) we can also detect the signal for 
smaller mass values of the doubly charged Higgs boson. 

Since we simulate signal and background events we define the signal significance as SS = 
Ns/y/Ng where Ns and Nb are the signal and background events in the chosen mass bin, 
respectively. Taking the corresponding integrated luminosities for ILC and CLIC e~e~ 
options L int = 10 4 pb _1 and taking into account the detection efficiencies for the leptons 
we calculate the number of events for the signal and background leading to the significances 
given in Table HVl Here, we find 5mu >> T H — , and therefore the signal significance is 
affected by the ISR+BS, FSR and smearing. If we consider lepton flavour conservation 
we can calculate the signal significance depending on the mass and coupling h ee of the 
doubly charged Higgs boson by applying appropriate px cuts. We may define the bins 
of width Am around the signal peak in the invariant mass distributions for each value of 
m H — . We then look at the fluctuations in the SM background in that bin and compare 
the corresponding rate for the signal. This procedure can be repeated for different couplings 
to find limiting values. The reaches of the coupling constant for the mass values of doubly 
charged Higgs boson smaller than the center of mass energies are given in Table IIV1 The 
resonance point worth a detailed study since smearing, initial and final state radiation affect 
the coupling dependence of the cross section. We estimate the reach of the coupling at which 
the resonances in the invariant mass distribution can be identified over the fluctuations in the 
SM background at 3<r level. We see from Fig. [HI that we can probe the doubly charged Higgs 
boson (when produced at resonance) coupling down to h ee = 3 x 10~ 4 and h ee = 10~ 3 with 
95% C.L. at the center of mass energy E cm = 500 GeV and E cm = 3000 GeV, respectively. 

For the process e~e~ — > ix~\iT the differential cross sections can be obtained as shown 
in Fig. HOI As can be seen from this figure two-muon invariant mass distribution shows a 
sharper peak than that of two-electron. In Fig. [TT] invariant mass distributions have been 
shown for different couplings at the resonance \fs = 500 GeV. In case of the lepton number 
violating process e~e~ — > ix~\f~ there is no significant SM background. We obtain the 95% 
C.L. contour by taking Poisson variable as the observed events with the mean of 9. In Fig. 
H2]the discovery contour of hu and m H — are plotted for the lepton number violating process 
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Figure 7: The invariant mass distribution of two electrons in the final state at the ILC energy 
yGJ = 500 GeV; (a-d) the upper four plots present the signal with m H — =350, 400, 450, 500 GeV. 
In (d) we show the signal distribution for the coupling h ee = 0.1. (e) the lower plot is for the 
background (ee — ► ee) with px e > 100 GeV and px e > 200 GeV. 

including ISR+FSR and smearing effects. As can be seen from Fig. [12] the accessible region 
of the contour plot can be enriched for higher integrated luminosities at future linear collider 
options. 

If we use an electromagnetic calorimeter with a larger energy resolution, we obtain a 
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Figure 8: The invariant mass distribution of two electrons in the final state at the CLIC energy 
^/i = 3000 GeV; (a-d) the upper four plots present the signal with m H — =500, 1000, 2000, 3000 
GeV. In (d) we show the signal distribution for the coupling h ee = 0.1. (e) the lower plot is for the 
background (ee — ► ee) with px e > 200 GeV and pr e > 500 GeV. 



smaller cross section for the signal and a larger one for the background (e~e~) to satisfy 
95% CL. Therefore, our results for S/yB will be changed accordingly and discovery region 
for doubly charged Higgs boson will become more restricted. 
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Table IV: Signal significances [SS(h ee ) = S(h ee )/yB] for the process e e — ► e e at the center of 
mass energy y^ = 500 GeV and y^ = 3000 GeV taking the integrated luminosity of L int = 10 pb -1 . 
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Figure 9: Signal significances depending on the coupling h ee for doubly charged Higgs boson (when 
produced at resonance) at the collider energies E cm = 500 GeV and 3000 GeV. 

V. CONCLUSIONS 

We have investigated the potential of the e~e~ running mode of the linear colliders for 
discovering and studying doubly charged Higgs bosons. In e~e" collisions the beams will 
deflect each other leading to a reduction in luminosity comparing to the e + e~ collisions in 
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Figure 10: Invariant mass distribution of the two-muon in the final state for different mass values 
of doubly charged Higgs boson at the center of mass energy t/s = 500 GeV. 

which two beams focus each other. For the part of the luminosity spectrum close to the 
nominal center-of-mass energy, the reduction in e~e~ mode is thus a factor of about 2.5 
compared with e + e~ |2J. The particularly clean environment of the e~e~ collisions justifies 
here the neglect of systematic errors. For luminosity higher than what we have used in our 
calculation the reach for couplings can be further enhanced. 

Using highly polarized electrons in the initial state make the resonances for doubly 
charged Higgs bosons observable down to quite smaller couplings. It has the benefit that at 
the same time one may also switch off most of the SM background processes (including W~ 
bosons) which could mask the discovery of new particles. The linear collider sensitivity to 
doubly charged Higgs boson properties is best when running at the resonance \/s = m H — . 
At the ILC and CLIC energies the couplings to charged leptons can be tested down to 
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Figure 11: Invariant mass distribution of two-muon around the resonance. The curves are obtained 
for different couplings hu by taking into account ISR+FSR and smearing effects at y/s = 500 GeV. 

3 x 10~ 4 and 10~ 3 at the resonance, respectively. Our analysis depends on fixed center of 
mass energy of the machine. Around the resonance there is still radiative return to the 
doubly charged Higgs boson resonance due to the initial state radiation, which makes it 
possible to discover these particles for small couplings (down to 5 x 10~ 2 for m H — = 1500 
GeV at \fs = 3000 GeV) without having a collider energy scan. The numerical analysis we 
performed is also relevant for any model which contains a Higgs triplet (e.g. Higgs triplet 
model, L-R symmetric models and little Higgs models). 

It is apparent that if a doubly charged Higgs boson is found at the LHC, an e~e~ collider 
would separately and in combination provide enormously important information concern- 
ing the structure and interactions of the Higgs sector. This also illustrates an important 
complementarity between the ILC or CLIC and the LHC. The discovery of a H~ " prior 
to the construction and operation of the e + e~/e~e~ collider would be very important in 
determining the range over which good luminosity and good energy resolution for e~e~ col- 
lisions should be a priority. Our calculations show that the signal significance SS \/Z can 
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Figure 12: Contour plot of the doubly charged Higgs boson coupling hu and mass m# — for the 
lepton flavor violating process e~e~ —* \x~~~ fj," at ILC and CLIC energies. 

be rescaled for a lower luminosity in case of a limited time of e-e- running mode of ILC and 
CLIC. 

On discovering a H it would be important to measure the absolute value of hij, there- 
fore a detailed experimental simulation of the decay modes is encouraged at both hadron 
and lepton colliders. 
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